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Abstract 

We solve the Ginzburg-Landau equation (GLE) for the mesoscopic super- 
conducting thin film of the square shape in the magnetic field for the wide 
range of the Ginzburg-Landau parameter 0.05 < Kg// < oo. We focus on the 
region of the field where formation of the antivortex has been reported previ- 
ously. We found that the phase with the antivortex exists in the broad range 
of parameters. When the coherence length decreases the topological phase 
transition to the phase with the same total vorticity and a reduced symmetry 
takes place. The giant vortex with the vorticity m = 3 is found to be unstable 
for any field, ^/a and /tgjj > 0.1. Reduction of Kgjj does not make the phase 
with antivortex more stable contrary to the case of the cylindric sample of 

the type I superconductor. 
PACS: 74.60.Ec 74.25.Ha 74.80.-g 
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I. INTRODUCTION 



Recently it was shown that the influence of boundaries can lead to stabilization of the 
vortex-antivortex molecules in mesoscopic samples Analysis of the linearized Ginzburg- 
Landau equation (GLE) has shown that such molecules appear at particular values of the 
external magnetic fleld depending on the sample shape and size 0. The solution of the 
GLE in the limit of the extreme type-II superconductor shows that such molecules have a 
very shallow minimum in the free energy and are very sensitive to the change of the 
sample shape 

In a square mesoscopic thin film with the total vorticity m = 3 the symmetric solution 
with four vortices and one antivortex is the solution of the linearized GLE with the lowest 
free energy According to ref. 0, away from the ifc2 line the giant vortex with vorticity 
m = 3 is stable and has the lowest free energy. This implies that a topological phase 
transition without change of the vorticity and without a reduction of the symmetry should 
take place with the change of the external field or /and the coherence length away from the 
critical-field line. 

It was proposed that in the limit k ~ l/\/2, where the vortex-antivortex interaction 
changes the sign, vortex-antivortex complexes should be more stable for the cylindric 
sample shape. For a thin film with Ag// = X^/d and Heff = ^eff/^, where d is the thickness 
of the film, A is the London penetration depth, and ^ is the superconducting coherence length, 
this mechanism should be less effective due to a smaller contribution of the 'magnetic energy' 
to the total free energy. 

We performed an extensive study of the region of the phase diagram where the vortex- 
antivortex phase was previously reported for the sample of the square shape. We focused 
to the region 4 < a/^ < 8 and Heff > 0.05. We found that the antivortex phase is stable in a 
broad range of parameters. The region of stability of the phase does not depend strongly on 
the value of the parameter Heff- The energy gain due to the antivortex formation is much 
smaller then the energy difference between two phases with different vorticities. The giant 
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vortex with m = 3 is unstable for any field, ^/a and k^jj ^0.1. Phase transition to the 
phase with three separated vortices takes place when ^/a is driven away from the critical 
field line. The reduction of Kg// does not stabilize the antivortex phase for the thin film 
sample in the contrast to the case of the cylindric sample ref. 

II. FORMALISM AND SOLUTION 



GLE for the normalized complex order parameter ■?/; = ^/^q, = \Jj3/\a\ has the 
following form: 

ei^y+'^-^f^-i^+m' = ^ (1) 

here ^ = 4^^^, o. and P are the temperature dependent parameters of the Ginzburg-Landau 
expansion for the free energy, $0 is the flux quantum, A is the vector potential and H = 
V X A the magnetic field. The second GLE equation for the vector potential reads: 



V X V X A = -^^(^*V^ - ^V^*) - (2) 



In addition to Eq.(l) we assume the boundary condition for the superconductor-insulator 
junction on the sample edges: 

(^V + ^) ■ n^; = 0, (3) 

where n is the vector normal to the surface of the sample. 

As it was described in ref. ^ we introduce N x N discrete points on the square and 
rewrite Eq.(l) in the form of the nonlinear discrete Schrodinger equation: 

Yl - e(i)ii,iV'i + V^ilV'iP = 0, (4) 

1 

where the summation index I = (±1,0), (0, ±1) points toward the nearest neighbours and 
^ii,i = i^N/ 0'Y^xp{'>'(Ph,i) (pi^^i = — ^ J^.'^ A(r)(ir. The boundary conditions are included 
in the discrete nonlinear Schr"odinger equation as in ref. where ■j/'i = if i is outside of 
the sample and e(i) = 4 - - 5i^^N - ^t^^i - dt^^N where \ = {i^ = I, . . . ,N,iy = I, . . . ,N). 
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After discretization of Eq. (2) we can obtain the exact expression for the vector potentiah 

A\ = Y.K{\-n)Jl (5) 

n 

where 

ea;p(-i0i-i,,i)V^i*V^i_i,) (6) 

where v G {x, y} and Ix = (1, 0), ly = (0, 1) 

Kin) = dxdy^^±d^±^ (7) 

27r^« -^0 _ 2cos (x) - 2cos (y) 

The numerical self consistent solution of the problem is obtained by iterating the solution 
of the nonlinear equation for the order parameter Eq.(4) and calculations of the current and 
the vector potential Eqs.(5,6). We used two ways of solving Eq.(4). The first is similar to 
that reported in ref. and corresponds to the iterative solution of the linearized Eq.(4). 
The second relies on the fact that Eq.(4) represents the Euler equation for the free-energy 
functional whith included boundary conditions. Eq.(4) was therfore solved by the direct 
minimization of the corresponding functional using the conjugate-gradient method. Both 
techniques gave identical results. 



III. RESULTS 

The main goal of the paper is to investigate the phase diagram in the region 4.5 < $/$o < 
6.5 and < 60 where the solution with one antivortex and four vortices (Fig.l) has 

been reported. We found that the region of the phase diagram where the symmetry induced 
antivortex solution has the lowest energy is broader than expected from the solution of the 
linearized GLE. As it is shown in Fig. 2 for n^jf = oo the antivortex phase is stable up to 
{(^/O"^ ~ 55, depending on $/$o- For a finite Keff this region shifts to the higher field as 
(a/^)^ increases (see Fig.2). 
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The interesting behaviour is observed when the external field is fixed and {a/^Y increases. 
Close to the Hc2 the lowest minimum of the free energy corresponds to the solution with the 
vorticity m = 4 — 1 with the antivortex in the center of the square. Present calculations do 
not confirm the existence of the giant-vortex solution with m = 3 in this region of the phase 
diagram as reported previously 0. The difference is due to increase of the number of discrete 
points enabling detection of the antivortex. With increase of {a/^Y away from the Hc2 
line the phase transition to the multivortex state with the same vorticity (m = 3) and a 
lower symmetry takes place (see Fig. 2). In general, the free energy depends on the vorticity 
m = n^ — ri- and the total number of vortices in the system n = 72+ + The transition at 
{a/^Y ~ 55 and $/$o ~ 5.5 takes place at the constant vorticity m = 3 with the change of 
n from 5 to 3. The transition is therfore not only characterized by an order parameter, but 
also by the change of the number of vortices at the constant total vorticity m, suggesting 
that the transition is close to the first order. This statement is confirmed by the observation 
that above the transition point, {a/S,) > {(i/Ocrit, both solutions with m = 3 and m = 4 — 1 
coexists (see Fig. 4). Since near the transition the free-energy difference between the phases 
with the same vorticity m and different n is small it is difficult to determine the phase 
boundary between phases with m = 4 — 1 and m = 3 accurately. The transition could be 
easier observed by calculating the two component order parameter rj^ = J x\ip{x,y)\'^dxdy, 

— i y\'^{^i y)\'^dxdy shown in Fig.3. The transition point is given by the point where the 
order parameter goes to 0. At the same point the calculated sample magnetization changes 
slope as it is clearly seen from Fig. 4. 

Close to the Hc2 line the repulsion of vortices from the boundaries and attraction of the 
4 vortices to the antivortex stabilizes the phase with m = 4 — 1 and small vortex-vortex 
distances. At smaller value of ^ the repulsion from the boundaries decreases and one vortex 
annihilates with the antivortex. As a result, the repulsion between the remaining vortices 
increases leading to an increase of the order parameter with a further decrease of ^. 

Increasing the field up to $/$o = 10.6 leads to the stabilization of the phase with total 
vorticity m = 7. Near the Hc2 line similar to the phase with m = 4 — 1 the solution with 
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m = 8 — 1 is realized. When {a/^Y increases in a complete analogy to the case with m = 4 — 1 
the second order phase transition to the phase with m = 7 and a similar order parameter 
takes place. Here also, both solutions with m = 8 — 1 and m = 7 coexists above {a/^)crit 
indicating that the transition is close to the first order. We believe that the situation is 
quite general for the case of arbitrary m = 4/ — 1 for / = 1, 2, 3, .... 

At the end we would like to discuss the dependence of the stability of the antivortex 
phase at small k. According to the arguments of Ref. at small k the vortex- vortex 
interaction changes the sign making the antivortex phase more stable. As a result, the 
average distance between vortices in the middle of the square increases as well. In order to 
verify this conjecture for the thin film sample we plot in Fig. 5 the vortex-antivortex distance 
tq as a function of l/^e//. The distance decreases with the decreasing ^e//- For K^jf < 0.1 
the distance is smaller than the grid spacing a/N so we can not resolve separate vortices. 
We find that tq oc exp (— A/Ae//) with A ~ a. The situation is just opposite to that reported 
in ref. [^]. We believe that in the case of the thin film of the square shape the reduction of 
K does not stabilize the phase with the antivortex. 

It is interesting to note differences between samples of different shapes. For the cylindric 
shape the giant vortex phase with any vorticity is always stable close to the Hc2 line 
According to the Ref. for the mesoscopic triangle the giant vortex state with m = 2 
is metastable and the solution with the antivortex (m = 3 — 1) is stable. For the case of 
the square shape the giant vortex solution with m = 3 is never stable for Kg// ^ 0.1. For 
Keff < 0.1 the limited grid spatial resolution prevented us to distinguish the solution with 
the antivotrex from the possibly (meta) stable giant- vortex solution. 

Finally, let us discuss the possibility to detect the state with the antivortex experimen- 
tally. Calculation of the magnetic field in the sample shows that the magnetic field has a 
local minimum in the center of the sample also for the giant vortex solution with m = 2. 
The local minimum observed for the antivortex state with m = 4 — 1 is therfore not due to 
the antivortex formation (Fig. 6) but due to a particular distribution of the current in the 
sample. Therefore, imaging of the magnetic field distribution cannot provide an evidence 
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for the antivortex. The magnetic field for the multivortex solution with m — 3 has 3 well 
separated maxima that break the four fold rotational symmetry of the sample allowing a 
direct imaging of vortices. Since the antivortex state cannot be detected directly the obser- 
vation of a hysteresis in the vicinity of the transition line from the m — A — 1 antivortex 
state to the m — 3 multivortex state could suggest that the symmetric phase is indeed the 
phase with the antivortex. 
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Figure captions 

Fig. 1 (a) The magnitude of the superconductivity order parameter for (a/^)^ = 35, 
Keff = oo and $ = 5.9<l'o- (b) The central region where the vortices are located in an 
expanded scale. The position of the antivortex is indicated by the symbol 0. 

Fig. 2 The calculated phase diagram. Different phases are marked with icons schemat- 
ically indicating the vortex pattern where the full dot represents a vortex, the open dot 
reperesents an antivortex and the larger full dot represents a double vortex. The full sym- 
bols and continous lines represent the phase boundaries for Keff = oo while the open symbols 
and dotted lines reperesent the phase boundaries for n^jf = 1. In the later case only the 
phase boundaries of the region with the total vorticity 3 are shown. 

Fig. 3 The magnitude of the order parameter around the transition where one vortex 
anihilates with the antivortex as a function of {a/^Y ^^e constant magnetic field. 

Fig. 4 The magnetic moment of the sample as a function of (a/0^ the constant 
magnetic field. 
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Fig. 5 The vortex-antivortex distance (a) and the magnetic moment of the sample (b) 
as functions of the parameter l/neff. In (a) error bars represent the grid spacing and the 
sohd hne the exponential fit discussed in the text. 

Fig. 6 The magnetic field in the film in the case of (a) giant vortex with m — 2, (b) the 
antivortex solution with m — A — 1 and (c) three separate vortices with m — 3. Here Hq is 
the external magnetic field. 
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